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A new stellar library at the near-IR spectral region developed for the empirical 
calibration of the Ca ii triplet and stellar population synthesis modeling is presented. 
The library covers the range AA 8348-9020 A at 1.5 A (FWHM) spectral resolution, 
and consists of 706 stars spanning a wide range in atmospheric parameters. We have 
defined a new set of near-IR indices, CaT*, CaT and PaT, which mostly overcome 
the limitations of previous definitions, the former being specially suited for the mea- 
surement of the Can triplet strength corrected for the contamination from Paschen 
lines. We also present a comparative study of the new and the previous calcium in- 
dices, as well as the corresponding transformations between the different systems. A 
thorough analysis of the sources of index errors and the procedure to calculate them is 
given. Finally, index and error measurements for the whole stellar library are provided 
together with the final spectra. 

Key words: stars: abundances - stars: fundamental parameters - globular clusters: 
general - galaxies: stellar content. 



1 INTRODUCTION 

With this paper we start a series dedicated to the empirical 
calibration of the near- infrared Call triplet. The ultimate 
aim of this work for us is to use the strength of the Ca lines 
in this spectral range to investigate the stellar content of 
early-type galaxies. However, we expect that researchers in 
different areas (e.g. starburst and active galaxies, globular 
clusters, and stellar astrophysics) can make use of some of 
the results and information that will be presented through- 
out these papers. These include a new stellar library, a set 
of homogeneous atmospheric parameters, empirical fitting 
functions describing the behaviour of the Ca triplet with 
the stellar parameters, and stellar population model predic- 
tions. 

Traditionally, elliptical galaxies have been thought to 
be a uniform class of objects, with global properties chang- 
ing smoothly with mass and hosting old and coeval stellar 
population. However, over the last decade, a growing body 
of evidence is indicating that the formation processes and 
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star formation histories of, at least, an important fraction 
of early-type galaxies is much more complex and hetero- 
geneous. The apparent age spread among elliptical galax- 
ies (Gonzalez 1993; Faber et al. 1995; j0rgensen 1999), the 
distribution of element abundances (Worthey 1998; Peletier 
1999; Trager et al. 2000a) and the interpretation of the scal- 
ing relations (like the colour-magnitude or Mg2-cr relations: 
Bower, Lucey & Ellis 1992; Bender, Burstein & Faber 1993: 
Pedraz et al. 1999; Terlevich et al. 1999; Kuntschner 2000: 
Trager et al. 2000b), are some of the main issues in the 
present debate about the evolutionary status of early-type 
galaxies. A major common obstacle to tackle with the pre- 
vious problems is how to disentangle age and metallicity 
effects in the integrated spectrum of a composite stellar pop- 
ulation. 

The measurement and interpretation of line-strength in- 
dices in the spectra of early-type galaxies has been a fun- 
damental tool to address the above topics, the Lick/IDS 
(Image Dissector Scanner) system (Worthey et al. 1994 and 
references therein) being the most widely used (but see also 
Rose 1994). The Lick system makes use of strong spec- 
tral features in the blue spectral range (from A 4100 A to 
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A 6300 A), although only the main absorption lines in the 
narrower AA 4800-5400 A range have been extensively mea- 
sured with high precision in the spectra of early-type galax- 
ies. More recently, new spectral indices have been defined 
and calibrated in the bluer region (Jones & Worthey 1995; 
Worthey & Ottaviani 1997; Vazdekis & Arimoto 1999; Gor- 
gas et al. 1999), with enlarged capability to decouple age 
and metallicity effects. However, the spectral range of in- 
terest is still quite narrow. It is obvious that enlarging the 
wavelength coverage will allow us to investigate the relative 
contributions of different stellar types to the composite spec- 
trum. As an example, a comparison between the mean ages 
derived using, as an age discriminant, features in different 
spectral regions, can provide important constraints about 
the star formation history of the galaxies (Sanchez-Blazquez 
et al. 2000). Although an important effort has been made 
to extend the stellar population studies to the ultraviolet 
region (e.g. Ponder et al. 1998), this is not the case for the 
near-infrared. In fact, it is remarkable that, although the 
features in the near-infrared spectral range where already 
included in the early analyses of the extragalactic old stellar 
populations (e.g. Spinrad & Taylor 1971), 30 years later, and 
despite the advent of modern CCD detectors, the potential 
of this spectral range to investigate the stellar population of 
early- type galaxies is still almost unexploited. 

The Call triplet is one of the most prominent features 
in the near-IR spectrum of cool stars and its potential to 
study the properties of stellar populations has been exten- 
sively acknowledged in the literature (see Section 2.2 for a 
review of previous works on the subject). However, a reliable 
analysis of the Ca triplet measurements in integrated spec- 
tra rests on the comparison of the data with the predictions 
of stellar population models. The accuracy of such predic- 
tions is highly dependent on the input calibration of the cal- 
cium line-strengths in terms of the main atmospheric stellar 
parameters (namely effective temperature, surface gravity 
and metallicity). The quality of this calibration has been 
the major drawback of previous stellar population models 
which have included predictions for the strength of the cal- 
cium triplet (Vazdekis et al. 1996; Idiart, Thevenin & de 
Freitas Pacheco 1997, hereafter ITD; Mayya 1997; Garci'a- 
Vargas, MoUa & Bressan 1998; Schiavon, Barbuy & Bruzual 
2000; Molla & Garcia- Vargas 2000). Previous calibrations 
have been either theoretical (based on model atmospheres), 
with their associated uncertainties, or based on empirical 
stellar libraries with a poor coverage of the atmospheric pa- 
rameter space. Apart from this major problem, there are 
other factors which have compromised, with different im- 
pacts, the reliability of these previous papers. They include 
stellar libraries with too few stars, problematic index def- 
initions (e.g., they cannot be used for all spectral types), 
uncertainties in the input stellar parameters (which trans- 
late into unknown uncertainties in the derived predictions), 
and illnesses in the fitting procedures (a proper statistical 
analysis of the computed coefficients is seldom followed). 
Throughout this work we will review and comment in detail 
all these problems and how we have coped with them. 

The main motivation of the present series of papers 
is, therefore, to provide a reliable calibration of the near- 
infrared Call triplet that makes it possible to accurately 
calculate the behaviour of the calcium line-strengths in stel- 
lar populations with a wide range of ages and metallicities. 



The outline of the series is as follows. In this first paper we 
introduce the new stellar library and define a new set of in- 
dices to quantify the strength of the Ca triplet. Paper II is 
dedicated to the determination of the input atmospheric pa- 
rameters for the stellar sample. These are some of the basic 
ingredients of the empirical fitting functions and the spec- 
tral synthesis analysis that will be presented in Paper HI 
and IV, respectively. In this latter paper, we will compute 
the predictions for single-burst stellar population models in 
two parallel ways: providing the Ca triplet strength using 
the fitting functions of Paper HI, and supplying spectral en- 
ergy distributions in the AA 8348-9020 A range (following 
Vazdekis 1999). 

In Section ^ of this paper we review previous work on 
the Call triplet behaviour and its applications to different 
astrophysical problems. Section ^ describes the new stellar 
library as well as an overview of the observation and reduc- 
tion procedures. A thorough discussion of the Call index 
definition is given in Section ^, where we define new im- 
proved Can indices which overcome the problems affecting 
previous definitions. This section also includes a compara- 
tive study of the sensitivity of the Ca ll indices to effects like 
the signal-to-noise ratio or the spectral resolution, and a set 
of calibrations for converting between different index defini- 
tions. In Section |^, an extensive analysis of index errors is 
presented, discussing the different sources of random errors 
and systematic effects. As a supplement of this section. Ap- 
pendix A provides accurate formulae for the computation 
of random errors in the index measurements, together with 
an analytical estimate of errors from signal-to-noise ratios. 
Finally, in Section ^ we refer the reader to a web page in 
which we provide a database containing the spectra for the 
whole library, an electronic table listing full information for 
each star, and public FORTRAN routines to compute the new 
indices and their associated errors. 



2 PREVIOUS WORKS ON THE NEAR-IR 
CA II TRIPLET 

In a pioneer work, Merrill (1934) (see also Wilson & Merrill 
1937) explored for the first time the near-infrared region of 
stellar spectra, identifying the most relevant absorption fea- 
tures and noting that the lines of Call AA 8498, 8542, 8662 
are prominent in types A-M and may have an interesting 
relationship to absolute magnitude. This tentative predic- 
tion was confirmed by the subsequent works of Keenan & 
Hynek (1945) and Parsons (1964). In particular, this last au- 
thor showed how the strength of the Ca lines (blended with 
Paschen lines for early-type stars) increases with absolute 
magnitude (i.e. being larger for supergiants) at fixed spectral 
type. Later works by O'Connell (1973), Anderson (1974), 
Cohen (1978), and others (see below) confirmed these find- 
ings, establishing the near-infrared Ca triplet as one of the 
most striking luminosity-sensitive feature in the spectra of 
F-M stars. 



2.1 Calibrations of the calcium line-strengths 

After these early studies, several authors have tried to model 
the behaviour of the Ca triplet with the basic atmospheric 
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parameters (Teff, logg and [Fe/H]). Two main lines of ap- 
proach have been foUowed, either using empirical stellar li- 
braries, or using the predictions of stellar atmosphere mod- 
els. 



2.1.1 Empirical calibrations 

Among the works following the empirical method, we must 
remark on that of Jones, Alloin & Jones (1984, hereafter 
J A J), in which they found that the Ca triplet strength 
strongly correlates with gravity (in the sense that it in- 
creases for giant and supergiant stars) . They also noted that 
the residuals of this relation correlate weakly with metallic- 
ity. This strong gravity dependence was independently con- 
firmed by Carter, Visvanathan & Pickles (1986) and Alloin 
& Bica (1989), although they did not agree in their conclu- 
sions about the metallicity dependence, being stronger in 
the latter study. 

An important step to quantify the g and [Fe/H] de- 
pendence was given by the comprehensive study of Di'az, 
Terlevich & Terlevich (1989, hereafter DTT). Using a stel- 
lar library with an enlarged range in gravity and metal- 
licity, these authors quantified the biparametric behaviour 
with logg and [Fe/H], noting that, in the high-metallicity 
range, the strength of the Ca triplet depends only on grav- 
ity, whilst, for low-metallicity stars, [Fe/H] is the main pa- 
rameter. Zhou (1991, hereafter ZHO) carried out a similar 
analysis, using a higher spectral resolution (2 A) and ex- 
tending the sample of stars by including cooler stars (up to 
M7). He found that the effect of temperature, found to be 
negligible in the previous work, was important at low Tcff. 
He also noted that the gravity dependence was stronger for 
giants than for dwarfs. 

Another relevant study is that of Malhk (1994, 1997). 
Through the analysis of an ample stellar library at high reso- 
lution (0.4 A), this author confirmed the strong dependence 
on logg (and showed that it was stronger as metallicity in- 
creases) and the milder dependence on [Fe/H] (but, remark- 
ably, more important for supergiant stars, in agreement with 
DTT). Unfortunately, he does not provide any fitting func- 
tion, restricting the potential application of these qualitative 
results. ITD have derived empirical functions which predict 
the strength of the Ca triplet in terms of the three atmo- 
spheric parameters. Using a library with a good coverage 
in metallicity but lacking a representative sample of super- 
giants, they conclude that the gravity dependence is not as 
strong as previously reported, and that metallicity is the 
main parameter. 

It is clear that there are some apparent contradictions 
among the conclusions of these papers. A critical analysis 
of these will be presented in Paper HI of this series, where 
we will show that the apparent discrepancies are mainly due 
to the different ranges of stellar parameters in the employed 
stell ar li braries, the diversity of index definitions (see Sec- 
tion 4.2), and differences in the fitting procedures (e.g. in 
most cases, a proper statistical approach was not applied). 



2.1.2 Theoretical calibrations 

The theoretical modelling of the Ca triplet lines is not sim- 
ple since the line cores are formed in the stellar chromo- 



sphere and, therefore, non-Local Thermodynamic Equilib- 
rium (NLTE) models are required. For this reason, the first 
attempts to analyze the sensitivity of the Ca strength to the 
atmospheric parameters were restricted to the line wings. 
In this sense. Smith & Drake (1987, 1990) concluded that 
the dependence on metallicity should be larger than that re- 
ported by the empirical studies (like J A J), and that effective 
temperatures had a non negligible effect. Erdelyi-Mendes & 
Barbuy (1991), also using LTE models, extended the analy- 
sis to cooler stars and included the significant contribution 
of molecular bands. They agreed with Smith & Drake in that 
metallicity, rather than gravity, was the main parameter, the 
[Fe/H] dependence being stronger for lower gravities. 

Using NLTE models, J0rgensen, Carlsson & Johnson 
(1992) carried out, for the first time, the computation of 
the line-cores, deriving full equivalent-widths (they showed, 
however, that the effects of departures from LTE are neg- 
ligible since the equivalent widths are dominated by the 
line wings). Concerning the sensitivity to stellar parameters, 
they found that the gravity dependence of the Ca triplet is 
infiuenced by T^s and [Fe/H] in a complicated way. They 
provided fitting functions for stars in the range 4000 K 
< Toff < 6600 K and -1.0 < [A/H] < +0.2 dex. 

Recently, Chmielewski (2000) has presented a compre- 
hensive analysis of the theoretical modelling of the Ca lines. 
He stresses that previous works did not take into account 
the contribution of hydrogen Paschen lines to the integrated 
equivalent widths, showing that this contribution is signif- 
icant for dwarfs hotter than 5800 K and giants with Tea 
above ~ 5500 K. Interestingly, this effect explains most of 
the discrepancies between the previous empirical and theo- 
retical works concerning the temperature dependence of the 
Ca lines. His main conclusions about the Ca line-strength 
sensitivity to the stellar parameters are that the gravity de- 
pendence only holds for logg < 3 dex, metallicity is a sig- 
nificant parameter in all cases, and temperature effects can 
only be neglected for dwarfs between 5000 and 6000 K. 

Unfortunately, these theoretical studies, being very im- 
portant to help to understand the behaviour of the Ca triplet 
and to check the empirical studies, are restricted to cool 
stars (typically between F and K) and, obviously, their con- 
clusions can not be extrapolated to later spectral types to 
construct stellar population synthesis models of, relatively 
old, stellar systems. 



2.2 Applications of the Ca triplet to the study of 
stellar systems 

Over the years, the near-IR Ca triplet has been used exten- 
sively to address a number of topics in different areas. 

2.2.1 Stellar astrophysics 

In the field of stellar astrophysics, the triplet, together with 
other spectral features in this spectral range, has been used 
for the spectral classification of stars (Sharpless 1956; Bouw 
1981; Kirkpatrick, Henry & McCarthy 1991; Ginestet et al. 
1994; Munari & Tomasella 1999). These features have also 
been widely used as chromospheric activity indicators (e.g. 
Linsky et al. 1979; Dempsey et al. 1993; Montes et al. 1998). 
Furthermore, its sensitivity to surface gravity has been ex- 
ploited to identify supergiants in the Galaxy (e.g. Garzon 
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et al. 1997) or in Local Group galaxies (Mantegazza 1992; 
Humphreys et al. 1988; Massey 1998). 

2.2.2 Globular clusters 

In spite of the first empirical papers that remarked that the 
Ca triplet was mainly gravity driven (see above), Arman- 
droff & Zinn (1988, hereafter A&Z), in a pioneering study, 
found that, in the integrated light of galactic globular clus- 
ters, the Ca equivalent width was well correlated with metal- 
licity, indicating that this feature could be a fair metallicity 
indicator for old, and approximately coeval, stellar popula- 
tions, provided that the metallicity was lower than solar. 

Subsequent works in this field concentrated in the study 
of individual red giants of globular clusters. Armandroff & 
Da Costa (1991, hereafter A&D) derived an empirical re- 
lation between the cluster metallicity and the "reduced" 
equivalent width W' = -E;W(CaT) -|-c (V - Vhb), where c is 
a constant term, EW (C&T) is the pseudo-equivalent width 
(typically the sum of the two strongest lines of the triplet), 
and Vhb is the V magnitude of the horizontal branch. Note 
that W' is a monotonic function of metallicity because grav- 
ity effects are removed by the last term of the equation. This 
relation has been recalibrated by several authors (e.g. Da 
Costa & Armandroff 1995; Geisler et al. 1995; Rutledge et 
al. 1997ab, hereafter RHS) and extensively used to derive 
metallicities of galactic globular clusters (Armandroff, Da 
Costa & Zinn 1992; Da Costa, Armandroff & Norris 1992; 
Buonanno et al. 1995; Suntzeff & Kraft 1996; Rosenberg 
et al. 1998), clusters and individual stars in the Magellanic 
Clouds (Olszewski et al. 1991; Suntzeff et al. 1992; Da Costa 
& Hatzidimitriou 1998; Cole, Smecker-Hane & Gallagher 
2000), and Local Group dwarf spheroidals (Suntzeff et al. 
1993; Smecker-Hane et al. 1999). A comprehensive analy- 
sis of the different methods to measure cluster metallicities 
using the the Ca triplet can be found in RHS. 

2.2.3 Active galaxies and extragalactic Hll regions 

One of the fields in which the study of the Ca triplet has 
had a wide application is that of active galaxies and star- 
burst regions. Terlevich, Di'az & Terlevich (1990a) studied 
the Ca triplet in a sample of normal and active galaxies, 
concluding that LINERS and Seyfert 2 objects exhibit Ca 
strengths equal or larger than those found in normal ellip- 
ticals (in contrast with the blue absorption lines, which are 
usually diluted in active galaxies). Following the conclusion 
of DTT that Ca equivalent widths above 9 A are only found 
in red supergiant stars, they interpret this result as a strong 
evidence for the occurrence of nuclear starbursts in the cen- 
tral regions of active galaxies. This same approach has been 
followed by a number of authors (Forbes, Boisson & Ward 
1992; Garcia- Vargas et al. 1993; Gonzalez Delgado & Perez 
1996ab; Heckman et al. 1997; Perez et al. 2000) to study the 
evolutionary status of the nuclei and circumnuclear regions 
of active galaxies through the presence of red supergiants. 

The Ca triplet has also been detected in extragalactic 
Hll regions (Terlevich et al. 1990b, 1996), being interpreted 
as due to the recent star formation burst, and it has been 
used to study the star formation history of extragalactic 
super-star clusters (Prada, Greve & McKeith 1994; Garci'a- 



Vargas et al. 1997; Gonzalez Delgado et al. 1997; Gallagher 
& Smith 1999; Goudfrooij et al. 2001). 

2.2.4 Early-type galaxies 

The first papers analyzing the Ca triplet in early-type galax- 
ies tried to make use of its gravity sensitivity to constrain the 
dwarf/giant ratio in the integrated light of old stellar popu- 
lations. The main debate at that time was focussed on the 
possible dwarf enrichment in the nucleus of M 31, as com- 
pared to its bulge or to low luminosity ellipticals like M 32. 
While Faber & French (1980) and, to some degree. Carter et 
al. (1986) favored a dwarf-enriched stellar population, Cohen 
(1978) and AUoin & Bica (1989) found that metallicity ef- 
fects and contamination by a molecular band could account 
for the apparent enrichment, mainly based on measurements 
of the strength of the Nal feature at A 8190 A. 

One interesting result from these first studies was that 
it was found that the strength of the Ca triplet did not vary 
much among early-type galaxies (Cohen 1979). Note that 
although Cohen (1978) found that the nucleus of M 31 ex- 
hibited stronger Ca lines than M 32, Faber & French (1980) 
reported the opposite behaviour. Later, Bica & AUoin (1987, 
hereafter B&A) measured the Ca strength in the spectra of 
62 galactic nuclei (from E to Sc) , concluding that the equiva- 
lent widths were not related to galaxy types or luminosities. 
This was confirmed by the studies of Terlevich et al (1990a), 
in which they found a small spread in the Ca strengths of 
14 normal galaxy nuclei, and Houdashelt (1995), who noted 
that the Ca equivalent widths of a sample of 34 early-type 
galaxies did not vary significantly among galaxies of differ- 
ent color or absolute magnitude. 

The apparent uniformity of the Ca measurement in el- 
liptical galaxies is somewhat surprising since: i) It contrasts 
with the metallicity dependence found for the stellar samples 
and the reported behaviour in galactic globular clusters (see 
above) , ii) Previous stellar population models predict a high 
sensitivity of the Ca triplet to the metallicity of metal-rich 
old stellar populations (Garcia- Vargas et al. 1998; Schiavon 
et al. 2000), and iii) It has been argued that, in contrast to 
other a-elements, Ca abundance is not enhanced compared 
to Fe in bright elliptical galaxies (O'Connell 1976; Vazdekis 
et al. 1997; Worthey 1998; MoUa & Garcia- Vargas 2000), 
and variations of the Fe line-strengths among ellipticals, al- 
though difficult to detect, are not negligible (e.g. Gorgas, Ef- 
stathiou & Aragon-Salamanca 1990; Gonzalez 1993; Davies, 
Sadler & Peletier 1993; Kuntschner 2000). Therefore, more 
reliable predictions of the Ca triplet behaviour in the in- 
tegrated stellar population of high metallicity systems and 
new high-quality measurements are hardly needed to clarify 
this point. 

This problem is closely related to the radial behaviour 
of the Ca strengths within galaxies. Given the measured 
metallicity gradients usually found in elliptical galaxies (see 
Gonzalez & Gorgas 1996, and references therein), even us- 
ing Fe line-strengths, one should expect to find unambiguous 
negative gradients in Ca. However, Cohen (1979) and Boro- 
son & Thompson (1991) found negligible gradients; Carter 
et al (1986) and Peletier et al. (1999) measured either null 
or positive gradients, and Delisle & Hardy (1992, hereafter 
D&H) found negative gradients in many galaxies of their 
sample. 
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Another important difficulty arises when trying to com- 
pare the absolute values of the Ca triplet measured in early- 
type galaxies with the actual values predicted by stellar pop- 
ulation synthesis models (e.g. Peletier et al. 1999; MoUa & 
Garcia- Vargas 2000). The uncertainties of the theoretical 
or empirical Ca behaviour implemented in such models, to- 
gether with other problems such as observational uncertain- 
ties, corrections from velocity dispersion or transformations 
between different index definitions, make very difficult to 
extract any significant conclusion from such a comparison. 



3 THE NEAR-IR STELLAR LIBRARY 

3.1 Previous libraries in the near-IR spectral 
range 

Previous stellar libraries providingspectra in the Ca ll triplet 
spectral range are listed in Table y. Here we only include li- 
braries with a spectral resolution better than 10 A (FWHM). 
For a list of lower resolution libraries we refer the reader to 
Munari & Tomasella (1999). We have also excluded from this 
list some libraries which are too specific for the purposes of 
this work, like those of Fluks et al. (1994) (only M stars), 
Allen & Strom (1995) (stars from two open clusters), RHS 
(huge collection of red giants from galactic globular clus- 
ters), or Montes & Martin (1998) (high resolution library in 
a limited spectral range). It is clear from this table that no 
previous library provides simultaneously a broad coverage 
of stellar metallicities and effective temperatures. In fact, 
the only libraries including a good fraction of low metallic- 
ity stars are those of DTT and ITD, but the spectral type 
ranges spanned in these cases are rather limited, especially 
in the latter case (as we will see in Paper IV, the relative 
contribution of stars colder than K3 to the near-IR spectrum 
of an old stellar population is very important). On the other 
hand, libraries composed of stars of all spectral types do not 
attain a broad range in metallicity (e.g. Munari & Tomasella 
1999). The low frequency of supergiant stars is an additional 
problem of some libraries (JAJ; ITD). It must be noted that 
some of the apparent discrepancies among different authors 
who have modelled the behaviour of the CaT index with 
the atmospheric parameters are mainly due to this different 
coverage of the parameter space by the calibrating stars (see 
Paper III). 



3.2 The sample 

We have observed a new stellar library of 706 stars at the 
near- infrared spectral region (AA 8348-9020 A). It includes 
421 of the 424 stars with known atmospheric parameters 
of the Lick/IDS Library (Burstein et al. 1984; Faber et al. 
1985; Burstein, Faber & Gonzalez 1986; Gorgas et al. 1993; 
Worthey et al. 1994). This subsample spans a wide range 
in spectral types and luminosity classes. Most of them are 
field stars from the solar neighbourhood, but stars covering 
a wide range in age (from open clusters) and with differ- 
ent metallicities (from galactic globular clusters) are also in- 
cluded. In order to obtain a large sample of stars in common 
with other previous works devoted to the Can triplet, our 
stellar library was enlarged to include 105 of the 106 stars 



of DTT and 43 of the 55 stars of the sample of ITD. More- 
over, with the aim of filling gaps in the parameter space, we 
also included stars from several other compilations, mainly 
focusing on O, B and A types (37 hot stars from the sample 
of Andrillat, Jaschek & Jaschek 1995), late M types, metal 
poor, metal rich and chemically peculiar stars (40 stars from 
a list kindly provided by G. Worthey, private communica- 
tion, and 102 stars from Jones 1997). The presence of hot 
stars in the stellar library allows us to expand the predic- 
tions of our models to young stellar populations. Also, and 
as it will be shown in Paper IV, very cool stars are neces- 
sary to reproduce the integrated spectra of old stellar sys- 
tems at the near-infrared spectral region over a wide range 
of ages. Finally, since calcium is an Q-element, 29 additional 
stars with high or low Ca/Fe ratios from the catalogue by 
Thevenin (1998) were observed in order to analyze the Call 
index dependence on relative abundances (see Paper III). 

The new stellar library covers the following ranges in 
atmospheric parameters: Tcff from 2750 to 38400 K, logg 
from to 5.12 dex, and [Fe/H] from -3.45 to +0.60 dex. It 
is important to note that the stars of the stellar library do 
not homogeneously cover the parameter space. In fact, most 
of the stars (~ 83 per cent) have metallicities from —1.0 to 
+0.5 dex, and the widest gravity range is also within this 
metallicity interval. The final atmospheric parameters which 
have been adopted for each star and the method to derive 
them are discussed in detail in Paper II of this series. We 
refer the reader to that paper for an HR diagram of the 
whole sample. 

3.3 Observations and data reduction 

The spectra of the stellar library were obtained during a to- 
tal 21 nights in six observing runs from 1996 to 1997 using 
the JKT (Jacobus Kapteyn Telescope), INT (Isaac New- 
ton Telescope) and WHT (William Herschel Telescope) at 
the Roque de los Muchachos Observatory (La Palma, Spain) 
and the 3.5-m telescope at Calar Alto Observatory (Almerfa, 
Spain). It is important to note that most of the stars (~ 93 
per cent), including all the field stars and the brightest ones 
from open and globular clusters, were observed during the 
three runs at the JKT using the same instrumental configu- 
ration. This ensures a high homogeneity a mon g the data of 
these three runs (see the analysis in Section 5.3). The instru- 
mental configuration at the other telescopes was also chosen 
to be as similar as possible to that at the JKT, obtaining 
spectral resolutions (FWHM) for the six runs in the range 
from 1.07 A to 2.13 A. A full description of these and other 
instrumental details is given in Ta ble tA We must note that, 
as it will be explained in Section B^^Tspectra from runs 1, 
3 and 4 were broadened to the spectral resolution of run 2 
(1.5 A). Hence, and with the exception of the few stars ob- 
served in runs 5 and 6, the whole library has a common 
resolution. A'obs stands for the number of stars observed in 
each run. 

Typical exposure times varied from a few seconds for 
bright stars to 1800 s for the faintest cluster stars. These 
ensured typical values of SN{A) ~ 100 A~^ (signal-to-noise 
ratio per angstrom) for field and open cluster stars, and 
SN(K) > 15 for the globular cluster stars. In order to test 
the derived random errors of the final index measurements 



(see Section 5T), we performed multiple observations of a 
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Table 1. Medium-resolution stellar libraries in the Can triplet spectral region. 



RiGfsrGIlCG 


No. stars 


R,Gsolu.tion 

I JT VV mVl , r\. J 


Spectral 
types 


[Fe/H] 
range 


CoiumGiitrS 




Jones et al. (1984) 


62 


3 


B - M5 


-0.60 , + 


0.23 






Diaz et al. (1989) 


106 


3.5 


F5 - Ml 


-2.70 , + 


0.55 






Zhou (1991) 


144 


2 


F5 - M7 


-2.28 , + 


0.60 


2 stars with [Fc/H] < - 


-0.61 


Andrillat et al. (1995) 


76 


1.2 


05 - GO 


7 








Serote Roos et al. (1996) 


21 


1.25 


B3 - M5 


-0.15 , + 


0.39 


Only giants and suporj 


jiants 


Carquillat et al. (1997) 


36 


2 


F5 - M4 


7 








Idiart et al. (1997) 


55 


~ 2 


Al - K3 


-3.15 , + 


0.35 






Mallik (1997) 


146 


0.4 


F7 - M4 


-3.0 , + 


1.01 


1 stars with [Fc/H]< - 


-1.6 


Munari & Tomasella (1999) 


131 


0.43 


04 - M8 


-0.54 , + 


0.30 






Monies et al. (1999) 


130 


~ 0.7 


FO - M8 


-2.74 , + 


0.31 


3 stars with [Fe/H] < - 


0.5 


This work 


706 


1.5 


06 - M8 


-3.45 , + 


0.60 







Table 2. Observational configurations 



Run Dates Telescope Spectrograph Detector Dispersion AA Slit width FWHM A^obs 

(A/pix) (A) (arcsec) (A) 



1 


20-26 Sep 1996 


JKT 1.0 m 


RBS 


TEK#4 


0.85 


8331- 


-9200 


1.5 


1.18 


264 


2 


15-19 Jan 1997 


JKT 1.0 m 


RBS 


TEK#4 


0.85 


8331- 


-9200 


1.5 


1.50 


248 


3 


25-29 Jun 1997 


JKT 1.0 m 


RBS 


TEK#4 


0.85 


8331- 


-9200 


1.5 


1.07 


362 


4 


13 Aug 1997 


INT 2.5 m 


IDS 


TEK#3 


0.85 


8185- 


-9055 


1.5 


1.28 


33 


5 


4 Aug 1996 


WHT 4.2 m 


ISIS 


TEK#2 


0.79 


8222- 


-9031 


1.2 


2.13 


12 


6 


17-18 Nov 1996 


CAHA 3.5 m 


TWIN 


SITe#4d 


0.81 


8300- 


-9912 


2.1 


2.11 


15 



subsample of stars (15 - 20) within each particular run and 
in common with other diflFerent runs. Finally, to perform a 
reliable flux calibration, several (3 or 4) spectrophotomet- 
ric standards were observed along each night at different 
zenithal distances. 

The reduction of the data was performed with 
ReD^eI (Cardiel 1999), which allows a parallel treatment 
of data and error frames (see more details in Section ^ 
and, therefore, produces an associated error spectrum for 
each individual data spectrum. We carried out a standard 
reduction procedure for spectroscopic data: bias and dark 
subtraction, cosmic ray cleaning, flat-fielding, C-distortion 
correction, wavelength calibration, S-distortion correction, 
sky subtraction, atmospheric extinction correction, spec- 
trum extraction, and relative flux calibration. We did not 
attempt to obtain absolute fluxes since both the evolution- 
ary synthesis codes and the line-strength indices only require 
relative fluxes. Cluster stars were also corrected from inter- 
stellar reddening using the color excesses from Gorgas et al. 
(1993) and Worthey et al. (1994) and the averaged extinc- 
tion curve of Savage & Mathis (1979). 

In order to optimize the observing time during the JKT 
runs, we decided not to acquire flat-fleld and comparison arc 
frames for each individual exposure of a library star. Con- 
cerning the flat-field correction, after checking that small 
variations of the CCD temperature do not affect the flat- 
field structure and considering that it exclusively depends 
on the position of the telescope, we obtained a complete set 
of flat-field exposures by pointing at a grid of positions on 
the dome (with a resolution of 30° in azimuth and 15° in 



t http:/ /www. ucm.es/info/Astrof/reduceme/reduceme. html 



altitude) during the day. Later, using the alt-azimuth coordi- 
nates of each star, the closest flat-fleld frame over the whole 
set was finally used as its own fiat-fleld calibration frame. It 
is important to highlight the importance of an accurate flat- 
field correction when the fringe pattern becomes relevant, 
as it is usually the case in this spectral range. Depending 
on the properties of the CCD, the fringe effect introduced 
high frequency structures with an amplitude of up to 7 per 
cent of the true flux, but they virtually disappeared by us- 
ing our procedure. Flux differences between the normalized 
fringe patterns of different flat-field frames allows an estima- 
tion of the uncertainties introduced in the fringe correction 
procedure. These are always below 0.6 per cent. 

Concerning the wavelength calibration, we only ob- 
tained comparison arc frames for a previously selected sub- 
sample of stars covering all the spectral types and luminos- 
ity classes in each run. The selected spectra were wavelength 
calibrated with their own arc exposures taking into account 
their radial velocities, whereas the calibration of any other 
star was performed by a comparison with the most similar, 
already calibrated, reference spectrum. This working proce- 
dure is based on the expected constancy of the functional 
form of the wavelength calibration polynomial within a con- 
sidered observing run (which has been widely tested by com- 
paring the derived polynomials for each run). In this sense, 
the algorithm that we used is as follows: after applying a 
test x-shift (in pixels) to any previous wavelength calibra- 
tion polynomial, we obtained a new polynomial which was 
used to calibrate the spectrum. Next, the calibrated spec- 
trum was corrected from its own radial velocity and, finally, 
the spectrum was cross-correlated with a reference spectrun 
of similar spectral type and luminosity class, in order to de- 
rive the wavelength offset between both spectra. By repeat- 
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ing this procedure, it is possible to obtain the dependence of 
the wavelength offset as a function of the test x-shifts and, as 
a consequence, we derive the required x-shift corresponding 
to a null wavelength offset. Uncertainties in the wavelength 
calibration are estimated in Section ^. 

It is also important to account for the presence of 
telluric absorptions, mainly at the red end of the spec- 
tral range (the strongest lines are locate at AA 8952, 8972, 
8980, 8992 A). Fortunately, since these strong H2O lines 
do not affect the Call triplet region and the observations 
were performed under dry conditions, we did not correct for 
their contamination and it is possible that our final spectra 
present unremoved features at the red end. An illustrative 
example of the effect of telluric lines in the Can triplet re- 
gion is given in Stevenson (1994) and Chmielewski (2000). 



4 CA II TRIPLET INDICES DEFINITION 

Line-strength indices have been widely defined to obtain ob- 
jective measurements of any relevant spectral feature. For 
the Ca triplet, there exist previous index definitions from 
other authors which were optimized to measure the lines for 
a narrow range of spectral types (mainly G and K-types). 
These previous indices, however, are not appropriate for all 
the spectral types. Hot and cold stars show strong spec- 
tral features which were not taken into account for those 
index definitions and, as a consequence, most of the calcium 
measurements for these spectral types become unrealistic. 
Moreover, some of the previous definitions have not been op- 
timized for composite stellar systems. Some of them require 
spectra with high signal-to-noise ratio, becoming unpractical 
for faint objects. Besides, some indices show a strong depen- 
dence on spectral resolution, or velocity dispersion broaden- 
ing. In this sense, one should not forget that, in the study 
of stellar populations, to measure true equivalent widths is 
not as important as to obtain robust index measurements, 
specially in spectra with low signal-to-noise ratios or with a 
wide range of spectral resolutions. 

In order to cope with the above problems, we have de- 
cided to int roduce some new Can triplet index definitions. 

presents an overview of the strongest spec tral 



Section 4.1 



features in the range from 8350 to 9020 A. In Section \L2 
we analyze the behaviour and limitations of previous Call 
triplet index definitions for di ffere nt spectral types. Th e new 

whereas Section 



indices are defined in Section 4.3 



4.4 



served to study the dependence of the new and old index 
definitions en the S/N ratio, the combined effect of spectral 
resolution and velocity dispersion broadening, and the flux 
calibration. Finally, in Section 4.5 we present comparisons 



and calibrations between different Ca 11 index systems. 



4.1 The AA 8350-9020 A spectral region 

Prior to any index definition, we have explored how the 
main absorption features in our wavelength range change 
along the spectral type sequence (Fig. ^). Apart from sev- 
eral atomic lines of intermediate strength, as those of Fe I 
(AA 8514.1, 8674.8, 8688.6, 8824.2 A), Mg i (8806.8 A) 
and Ti I (8435.0 A), the Can triplet (AA 8498.02, 8542.09, 
8662.14 A) is the strongest feature over a wide range of 



spectral types (from F5 to M2 approximately) and for all 
luminosity classes (see Fig. ^c). 

The hydrogen Paschen series (AA 8359.0, 8374.4, 8392.4, 
8413.3, 8438.0, 8467.3, 8502.5, 8545.4, 8598.4, 8665.0, 
8750.5, 8862.8, 9014.9 A, from P22 to PIO respectively) is 
apparent in stars hotter than G3 types. Depending on the 
luminosity class, the strength of this series reaches a max- 
imum for F or A type stars (Andrillat et al. 1995). In the 
earliest spectral types, where the Ca 11 strength becomes in- 
significant, the relative depths of the Paschen lines show a 
smooth sequence with wavelength (see Fig. |l^). However, 
due to the fact that, for low and intermediate spectral res- 
olution, the Paschen lines P13, P15 and P16 overlap with 
the Can triplet, these three Paschen lines stand out in the 
smooth sequence for A and F types (Fig. |l|3). The fact that 
the Ca n lines are blended with the Paschen lines has always 
been an obstacle to measure the calcium triplet in warm 
stars (see B&A and Chmielewsky 2000). 

Stars cooler than early M types exhibit molecular bands 
that change the slope of the local continuum (Fig. ^) . The 
strongest ones are the triple-headed band at AA 8432, 8442, 
8452 A and the double-headed bands at AA 8859.6, 8868.5 A 
and AA 8937.4, 8949.8 A of titanium oxide (TiO). Other 
molecular features are the band sequence of TiO at AA 8373, 
8386, 8420, 8457, 8472, 8506, 8513, 8558, 8569 A and several 
vanadium oxide (VO) bands at A A 8521, 8538, 8574, 8597, 
8605, 8624, 8649, 8668 A (see Kirkpatrick et al. 1991 and 
references therein) . The strength of these features increases 
as the temperature decreases, being more prominent for gi- 
ants than for dwarfs. The spectra of late M type stars are 
dominated by strong molecular bands showing very weak 
Can lines. 



4.2 Previous index definitions 

Several previous works have established different index def- 
initions to measure the strength of the Can triplet. In par- 
ticular, the most commonly used indices are those defined 
by JAJ, B&A, A&Z, DTT, A&D, ZHO, D&H, and RHS. 
Most of these indices were defined according to the classical 
definition of line-strength indices, that is, by establishing a 
central bandpass covering the spectral feature of interest, 
and two other bandpasses (at the red and blue sides of the 
central region) which are used to trace a local continuum 
reference level through a linear fit to the mean values in 
both bands (for more details about the definition and com- 
putation of classical indices see Appendix Al). In Table ^ 
we list the bandpasses limits of the previous definitions. 

It must be noted that some of these indices slightly dif- 
fer from the classical definition. On the one hand, the con- 
tinuum by ZHO is computed from the mean value of the 5 
highest pixels in each continuum bandpass, whereas in the 
index by JAJ, it is chosen relative to the maximum fiux 
around two selected wavelengths. It is important to stress 
that this kind of index definitions is useless for low signal-to- 
noise spectra (S/N), leading to a spurious negative correla- 
tion between index value and S/N ratio. Also, since the clas- 
sical index definition refers to the measurement of a unique 
spectral feature, previous definitions are given for each Ca n 
line. In this sense, the most general Can triplet index is 
usually computed as the sum of the three single- line indices, 
although some authors (DTT, A&D and ZHO) prefer to use 
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Figure 1. Spectra of the stars HD 186568 (B8 III), HD 89025 (FO III), HD 216228 (KO III) and HD 114961 (M7 III) in the spectral 
range of the stellar library. The strongest features in this region are marked: the Paschen Series (from Pll to P20), the Call triplet (Cal, 
Ca2 and Ca3), several metallic lines (Fel, Mgl and Til), molecular bands (TiO and VO) and telluric absorptions. 



the sum of the two strongest Unes of the triplet (AA 8498.0 
and 8542.1 A) or even a weighted combination of the three 
hnes (RHS) in order to optimize the sensitivity of the index 
to the signal-to-noise ratio. 

Fig. ^ illustrates both the bandpass position and the 
predicted continuum of some definitions over three different 
spectral types. The indices by B&A and RHS have not been 
included since their blue bandpasses are out of our spectral 
range. Also the index by A&D has been excluded due to its 
similarity to the index by A&Z. 

In order to predict a reliable local continuum, the con- 
tinuum bandpasses should be located in spectral regions free 
from strong absorption features. In this sense, all the pre- 
vious definitions are optimized for G and K spectral types 
(see Figs. |2k and ^). However, because the rel- 

ative contribution of the earliest spectral types to the inte- 
grated spectrum of stellar systems was expected to be neg- 
ligible at this spectral range, most of the previous calcium 
indices were defined without taking into account the wave- 
length position of the Paschen series. The indices by A&Z, 
DTT, and D&H clearly suffer from such a limitation deriv- 
ing a continuum level below the true one for hot stars (see 
Figs. and ^m). On the other hand, because of the 

presence of molecular bands in late type stars, some defi- 
nitions show unreliable continuum levels for these spectral 
types. In particular, the red bandpass of the index by JAJ 
falls in a strong TiO absorption (Fig. ^), while the blue 



bandpasses by DTT and ZHO (which are roughly the same) 
are located too close to the continuum break, causing a rise 
in the continuum level at the position of the Ca lines (see 
Figs, ^i and |2|l) . Also, when the relative distance between the 
two continuum bandpasses is too large, the index can be in- 
sensitive to local variations in the true continuum (JAJ and 
DTT) , whereas indices with continuum sidebands which are 
too close (A&Z) are hig hly se nsitive to velocity dispersion 
broadening (see Section 4.4.2). As far as the width of the 
sidebands is concerned, bandpasses that are too wide (D&H) 
are not appropriate when the spectral region is full of ab- 
sorption lines which decrease the continuum level. Indices 
with two narrow sidebands (JAJ) are very sensitive to the 
S/N ratio and to residuals from sky subtraction. Another 
important factor is the width of the central bands, since, as 
we will see in Section 4.4, it mainly determines the sensitiv- 



ity of the index to the S /N ratio and the velocity dispersion 
broadening. Finally, and not the least important, it is ob- 
vious that any Ca index definition following the classical 
system is unavoidably affected by Paschen contamination in 
the central bandpass for stars hotter than Gl and G3, for 
dwarfs and giants respectively. 



4.3 New indices definitions 

Although the previous index definitions have been very use- 
ful to increase the understanding of the behaviour of the 
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Figure 2. Previous calcium index definitions over different spectral types. Indices (from top to bottom) correspond to the systems of 
JAJ, A&Z, DTT, ZHO and D&H. The representative spectra (from left to right) are those of HD161817 (A2 VI), HD25329 (Kl Vsb) 
and HD148783 (M6 III). Grey and open bands mark, respectively, continuum and central bandpasses, whereas the solid line represents 
the local pseudo-continuum computed by an error weighted least-squares fit to all the pixels in the continuum bands. 



Ca II triplet both in individual stars and stellar populations 
(see Section P), in this work we have decided to define a new 
set of improved indices, specially designed to be measured in 
the integrated spectra of galaxies. With the new definitions 
we have tried to alleviate difficulties such as the continuum 
definition, the Paschen contamination, and the sensitivity 
to the S/N ratio and the velocity dispersion (spectral res- 
olution). Note that the aim of these new indices is not to 



provide measurements as close as possible to the true equiv- 
alent widths. Other indices, like those of RHS are much more 
appropriate for that purpose. 

The new indices have been defined according to a new 
type of line-strength index concept: the generic index. It is 
a natural generalization of the classical definition which in- 
cludes the following requirements: it is characterized by an 
arbitrary number of continuum and spectral-feature band- 
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Table 3. Bandpasses limits of previous calcium indices. Codes are 
the following: JAJ (Jones, AUoin & Jones, 1984), B&A (Bica & 
Alloin, 1987), A&Z (Armandroff & Zinn, 1988), DTT (Di'az, Ter- 
levich & Terlevich, 1989), A&D (Armandroff & Da Costa, 1991), 
ZHO (Zhou 1991), D&H (Delisle & Hardy, 1992), and RHS (Rut- 
ledge et al., 1997a). Due to the subjective continuum definition 
by JAJ, bandpasses limits for this index are as given in D&H. 



Index 



Central 
bandpass (A) 



Continuum 
bandpasses (A) 



Cal(JAJ) 
Ca2(JAJ) 
Ca3(JAJ) 

Cal(B&A) 
Ca2(B&A) 
Ca3(B&A) 

Cal(A&Z) 
Ca2(A&Z) 
Ca3(A&Z) 

Cal(DTT) 
Ca2(DTT) 
Ca3(DTT) 

Cal(ZHO) 
Ca2(ZHO) 
Ca3(ZHO) 

Ca2(A&D) 
Ca3(A&D) 

Cal(D&H) 
Ca2(D&H) 
Ca3(D&H) 

Cal(RHS) 
Ca2(RHS) 
Ca3(RHS) 



8483.0-8511.0 
8517.0-8559.0 
8634.0-8683.0 

8476.0-8520.0 
8520.0-8564.0 
8640.0-8700.0 

8490.0-8506.0 
8532.0-8552.0 
8653.0-8671.0 

8483.0-8513.0 
8527.0-8557.0 
8647.0-8677.0 

8488.0-8508.0 
8532.0-8552.0 
8652.0-8672.0 

8532.0-8552.0 
8653.0-8671.0 

8483.0-8511.0 
8517.0-8559.0 
8634.0-8683.0 

8490.0-8506.0 
8532.0-8552.0 
8653.0-8671.0 



8633.0-8637.0, 8903.0-8907.0 
8633.0-8637.0, 8903.0-8907.0 
8633.0-8637.0, 8903.0-8907.0 

8040.0-8160.0, 8786.0-8844.0 
8040.0-8160.0, 8786.0-8844.0 
8040.0-8160.0, 8786.0-8844.0 

8474.0-8489.0, 8521.0-8531.0 
8521.0-8531.0, 8555.0-8595.0 
8626.0-8650.0, 8695.0-8725.0 

8447.5-8462.5, 8842.5-8857.5 
8447.5-8462.5, 8842.5-8857.5 
8447.5-8462.5, 8842.5-8857.5 

8447.0-8462.0, 8631.0-8644.0 
8447.0-8462.0, 8631.0-8644.0 
8447.0-8462.0, 8631.0-8644.0 

8474.0-8489.0, 8559.0-8595.0 
8626.0-8647.0, 8695.0-8754.0 

8559.0-8634.0, 8683.0-8758.0 
8559.0-8634.0, 8683.0-8758.0 
8559.0-8634.0, 8683.0-8758.0 

8346.0-8489.0, 8563.0-8642.0 
8346.0-8489.0, 8563.0-8642.0 
8563.0-8642.0, 8697.0-8754.0 



passes, the contribution of each spectral-feature bandpass 
can be modified by defining a multiplicative factor for each 
one, and the pseudo-continuum is derived by using an er- 
ror weighted least-squares fit to all the pixels of the con- 
tinuum bandpasses. Although at first sight it could seem 
that this approach is exactly the same that the one obtained 
by adding several classical indices (weighted with the corre- 
sponding multiplicative factors), the fact that all the spec- 
tral features share the same continuum has an important ef- 
fect in the computation of the index error as it is explained 
in Appendix A2. It is also convenient to highlight that the 
requirements incorporated in the definition of generic indices 
translate into important improvements: 

(i) Generic indices are specially suited for the measure- 
ment of adjacent spectral features, where the use of a com- 
mon continuum level can be useful. They are also highly 
recommended when the spectral region of interest is densely 
populated by spectral features, since many thin continuum 
bandpasses at adequate locations make it possible to avoid 
the presence of the other spectral features. Moreover, a large 
number of continuum bands guarantees a low sensitivity to 
the S/N ratio and a robust continuum definition. 

(ii) As we will prove later, the multiplicative factors are 
useful to remove the contamination by other spectral fea- 
tures (either in absorption or in emission). 



Table 4. Bandpasses limits for the generic indices CaT and PaT. 

CaT central PaT central Continuum 

bandpasses (A) bandpasses (A) bandpasses (A) 



Cal 8484.0- 


-8513.0 


Pal 8461.0- 


-8474.0 


8474.0- 


-8484.0 


Ca2 8522.0- 


-8562.0 


Pa2 8577.0- 


-8619.0 


8563.0- 


-8577.0 


Ca3 8642.0- 


-8682.0 


Pa3 8730.0- 


-8772.0 


8619.0- 


-8642.0 










8700.0- 


-8725.0 










8776.0- 


-8792.0 



(iii) The use of an error weighted least-squares fit in the 
determination of the local pseudo-continuum is specially ad- 
vantageous when we measure the near-IR absorption fea- 
tures. In this spectral region, the presence of sky emission 
lines and telluric absorptions implies that the signal-to-noise 
ratio, as a function of wavelength, is a highly inhomogeneous 
function. 

Following the notation from Cohen (1978) for the cal- 
cium triplet, the new generic index will be referred as CaT. 
It is defined by establishing five continuum bandpasses, 
and three central bandpasses covering each calcium line 
(Cal, Ca2 and Ca3). The continuum regions were carefully 
chosen to optimize the continuum level for all the spec- 
tral types, even when the spectra were broadened up to 
a = 300 km s~^. In order to calibrate the CaT contam- 
ination by the Paschen series in early spectral types, we 
defined a new generic index, namely PaT, which measures 
the strength of three Paschen lines free from the calcium 
contamination. The continuum bandpasses are the same as 
in CaT whereas the spectral- feature bandpasses (Pal, Pa2 
and PaS) are centered on the lines P17, P14 and P12 of the 
series. In both definitions, the multiplicative factors equal 
unity and the analytical expressions are 



CaT = Cal -f Ca2 -f Ca3, 
PaT = Pal -I- Pa2 4- PaS. 



and 



(1) 
(2) 



The bandpasses limits for the CaT and PaT indices are 
listed in Table W, and the bandpasses positions and predicted 
continua are illustrated in Fig. ^ for different spectral types. 

It is important to note that, apart from rough correc- 
tions applied to the integrated spectra of young star clusters 
(Alloin & Bica 1989) and theoretical considerations on syn- 
thetic spectra (Chmielewski 2000), few previous works have 
faced the Paschen contamination in individual stars. 

Finally, we have defined a new calcium triplet index, 
namely CaT* , which expresses the strength of the Ca ll lines 
corrected from the contamination by Paschen lines. The new 
index was defined by imposing the following requirements. 
First, the values of CaT* should be very similar to those 
of CaT for late-type stars. Secondly, since the true calcium 
strength in hot stars is very low or even null, the CaT* index 
should tend to values around zero for the earliest spectral 
types. Note that, since the CaT index in hot stars is actu- 
ally measuring the strength of the Peischen lines rather than 
those of Ca, such an index (as well as those used in previous 
works) can lead to wrong conclusions when interpreting the 
integrated spectra of young stellar populations. 

To estimate the hydrogen contribution to the calcium 
index, we have compared the CaT and PaT indices for hot 
stars, which show a pure smooth Paschen series in their spec- 



© RAS, MNRAS 000, [^-^ 



The near-IR Call triplet: stellar library 11 





8400 8500 8600 8700 

wavelength (A) 



Figure 3. CaT and PaT indices over different spectral types. The 
spectra correspond (from top to bottom) to HD161817 (A2 VI), 
HD25329 (Kl Vsb) and HD148783 (M6 III), which are the same 
spectral types displayed in Fig. ^. Grey and open bands represent 
continuum and central bandpasses respectively, with the solid line 
showing the derived local pseudo-continuum. 



tra. For this subsample, we should expect a hnear relation- 
ship between the CaT index (which measures the P16, P15 
and P13 lines falling in the Cal, Ca2 and Ca3 bandpasses) 
and the PaT index (measuring the P17, P14 and P12 lines). 
Such a relation was calibrated by using an error weighted 
least-squares fit to the 26 library stars with T^s > 10500 K 
(see Fig. |^ obtaining 

CaT = (0.93 ± 0.03) PaT. (3) 
Thus the Paschen corrected CaT* index was defined by sub- 




PaT(A) 

Figure 4. CaT - PaT diagram for the whole stellar library. Dif- 
ferent symbols are used to indicate different ranges of effective 
temperature, as shown in the key. The solid line CaT = 0.93 PaT 
represents an error weighted least-squares fit to the 26 filled sym- 
bols. 

tracting from CaT the correction term given by the above 
equation, i.e. 

CaT* = CaT - 0.93 PaT, (4) 

which can be considered, in turn, as a generic index with five 
continuum bands, six central bandpasses and two different 
multiplicative factors (1.00 and —0.93). Since we are adopt- 
ing a definition, the multiplicative factor —0.93 is assumed 
to be exact, so we are not concerned with the nominal error 
0.03 obtained in the fit to the hot stars. 

Finally, it is worth noting that CaT* is also useful for 
late type stars where PaT takes values around zero and the 
correction term becomes negligible. 

For those readers interested in measuring the new in- 
dices in their spectra, we provide a set of public FORTRAN 
subroutines available from the website given at the end of 
the paper. 

4.4 Sensitivities of the indices to different effects 

In this section, we discuss the sensitivity of the new and 
previous Ca indices to the signal-to-noise ratio, velocity dis- 
persion broadening (or spectral resolution), relative flux cal- 
ibration and sky subtraction. 

4-4- i Signal-to-noise ratio 

In Appendix A we show in detail how to compute random 
errors, arising from photon statistics, for the new generic 
indices. Following a procedure similar to that presented in 
Cardiel et al. (1998, hereafter CGCG), we have derived ac- 
curate formulae for the computation of random errors for 
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Table 5. Ratio between the relative random error of any index 
and that of CaT. 



(T[CaT]/CaT 



CaT(ZHO) 


0.788 


± 0.004 


CaT{D&H) 


1.009 


± 0.006 


CaT(JAJ) 


3.938 


± 0.046 


CaT(DTT) 


0.933 


±0.007 


CaT(A&Z) 


0.567 


±0.002 



the new Ca indices (Appendix A2). Note that reliable ran- 
dom errors can only be derived when, after a full control 
of the error propagation throughout the data reduction, an 
accurate error spectrum is derived for each data spectrum. 
Appendix A3 also provides analytical estimates of the pre- 
dicted random errors as a function of signal-to- noise (S/N) 
ratio per angstrom of the form 

Cl — C2 Xa 



SN{k) 



(5) 



The exact values of the Ci coefficients are given in 



Eqs. (A37). As an illustration, for an MO giant, a S/N per 
angstrom of 16 is required to measure the CaT and CaT* 
indices with a 10 per cent uncertainty. We have selected this 
spectral type since, as we will see in Paper IV, old stellar 
populations exhibit near-IR spectra close to those or early- 
M stars. 

To compare the sensitivity of our indices to the S/N 
ratio with that of previous definitions, we have computed 
the typical random error due to photon noise (o"[Ia]) for 
the 706 spectra of the library and all the considered index 
definitions. Table ^ lists the ratio between the relative ran- 
dom error of any index and that of CaT. The relations were 
computed through a least-squares linear fit (rejecting data 
outside the 99.73 per cent confidence level) to the relative 
errors of the whole sample. It can be seen that our index is 
similar to that of D&H, but leads to somewhat larger ran- 
dom errors than the definitions of ZHO, A&Z and DTT. 
The reason for this is that these systems use relative narrow 
central bandpasses (with a total width for the Ca2 and Ca3 
lines of 38, 40 and 60 A respectively, compared to the 80 A 
in our definition). Note that since the noise in the bottom 
of any absorption line is lower than in the wings, indices 
with narrow central bandpasses yield small random errors. 
Unfortunately, this turns into a high sensitivity to veloc- 
ity dispersion broadening (as we will see in Section 4.4.2). 
On the other hand, the index from JAJ uses a wide central 
band and extremely narrow continuum bandpasses which 
increase the uncertainties in the determination of the con- 
tinuum level. This is why this index shows a high sensitivity 
to the signal-to-noise ratio. 

As it was also noted in Section some Ca indices 
have been defined as the sum of the two strongest Ca ll lines 
in order to increase the S/N ratio of the final measurement. 
In our case, when the Cal line is removed from the CaT def- 
inition, the relative random errors decrease by a factor of 10 
per cent. However, keeping in mind that the central band- 
pass of the Ca3 line falls near the region of strong telluric 
absorptions, being completely immersed in it for galaxies 
with z > 0.03, we have decided to include the bluest line 



of the triplet (Cal) in the index definition, at the expense 
of sacrificing the gain in accuracy, to ensure a more reliable 
index measurement for galaxies with relatively high radial 
velocities. 



4-4-2 Spectral resolution and velocity dispersion 
broadening 

In order to study the sensitivity of the Ca indices to the 
spectral resolution or velocity dispersion (a), we selected 
a subsample of near-solar metallicity templates of all the 
spectral types and luminosity classes and broadened their 
spectra with additional a's from a — 25 km s^^ up to 
a — 400 km s~^ (in steps of 25 km s~^). The indices were 
measured in the full set of broadened spectra, and we com- 
puted, for each group of templates with the same spectral 
type, a third-order polynomial fit to the relative changes of 
the index values as a function of velocity dispersion 

la(o-) - Xa(o-o) , , , 2 , , 3 _ , 



(6) 



where ctq = 22.2 km s~^, corresponding to the final spec- 
tral resolution of the stellar library (FWHM — 1.50 A, see 
Section ^^). For a given spectral type, we did not com- 
pute separate fits for the different luminosity classes since 
the derived polynomials were very similar. Using the above 
equation, it is easy to show that the indices measured at two 
different velocity dispersions, cri and cr2, are related by 



Xa(cr2) = I'a(o-l) 



l-p(ai) 



(7) 



1 -p(o-2) ' 

Table ^ lists the derived coefficients for CaT* , CaT and 
PaT measured over different spectral types. Note that we 
only provide calibrations for spectral types with equivalent 
widths significantly different from zero (e.g. stars cooler and 
hotter than mid-F types for CaT* and PaT respectively). 
For illustration, such polynomials are shown in Fig. ^a. It is 
readily observed that the sensitivity of CaT* to the veloc- 
ity dispersion broadening is much lower than that of CaT 
and PaT. Such an improvement is explained since the effects 
of the broadening on CaT and PaT are partially compen- 
sated when the CaT* index is computed. This fact should 
be considered as another additional advantage of employing 
generic indices for the measurement of spectral features. 

The measurements and fitting procedures given above 
were also performed for all the previous Call indices, al- 
though only the c alibrations computed for MO type spectra 
(see Section 4.4.1) are presented in this paper. The derived 
coefficients are presented in Table ^, whilst the polynomials, 
together with those of CaT and CaT* are shown in Fig. ^d. 
It is clear that, except for the index by D&H, all the previ- 
ous indices are more sensitive to broadening than our defi- 
nitions. In most cases, the high sensitivity is due to the rel- 
atively narrow central bandpasses (A&Z, ZHO and DTT). 
The index by JAJ has wide central bandpasses, but the blue 
continuum band is located too close to the Ca3 line. Finally, 
the index by D&H, with wide continuum and central band- 
passes, is the only one which overcomes the sensitivity to 
broadening of CaT, with a similar behaviour to our CaT* 
index. 
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Figure 5. Index sensitivity to the spectral resolution or velocity dispersion broadening. Fig. ^ shows the relative variation of CaT* 
(solid lines), CaT (dashed lines) and PaT (dotted lines) for several broadened spectral types. In Fig. similar polynomials are shown 
for the new and the previous Call indices for an MO type star. In both diagrams, A/// is zero for a = 22.2 km s~^ (see Eq. (^)), 
corresponding to the spectral resolution of the stellar library (FWHM = 1.50 A). 



Table 6. Coefficients of the broadening correction polynomials 
(AXa/Xa = a + b a + c a'^ + d a'"^) ior the new indices and 
different spectral types. 



Xa SpT a(xl0-3) 6(xlO-^) c(xl0-7) d(xl0-9) 



CaT* F7 


0.784 


-4.683 


6.197 


-4.563 


G3 


-0.452 


1.G48 


5.386 


-4.201 


KG 


-1.565 


5.605 


7.483 


-4.418 


MO 


-6.24G 


27.524 


3.468 


-3.797 


M5 


-6.795 


27.98G 


12.715 


-3.986 


CaT AO 


5.562 


-23.662 


-4.621 


-7.408 


F7 


0.503 


-2.119 


G.559 


-5.523 


G3 


-1.969 


10.305 


-5.495 


-4.407 


KG 


-3.106 


15.413 


-5.453 


-4.332 


MO 


-8.542 


41.116 


-11.118 


-3.485 


M5 


-3.953 


17.938 


G.243 


-3.751 


PaT AO 


3.943 


-16.146 


-6.162 


-4.991 


F7 


2.169 


-7.226 


-10.861 


-2.664 


Table 7. Coefficients of the broadening correction polynomials 


(AXa/Xa = a 


+ b 0- + c 0-2 


+ d a^) for previous index definitions 


and an MG spectral type. 








Xa 


a(xl0-3) 


fe(xlG-5) 


c(xlG~^) 


d(xlG-9) 


CaT(A&Z) 


-7.712 


48.585 


-59.578 


-12.551 


CaT(D&H) 


-3.517 


15.486 


2.396 


-3.610 


CaT(DTT) 


7.1G7 


-31.378 


-2.G28 


-3.755 


CaT(JAJ) 


-49.068 


237.3G4 


-74.233 


4.143 


CaT(ZHO) 


-15.283 


83.299 


-65.193 


G.314 


4.4.3 Flux 


calibration 








There is no simple recipe to estimate the sensitivity of an in- 



dex to the uncertainties in the relative flux calibration, since 
it not only depends on the index deflnition but also on the 
response curve of the instrumental configuration. In general, 



the effect is larger as the wavelength coverage of the side- 
bands increases. As an illustration, for the JKT runs, the 
systematic differences between the CaT index measured af- 
ter and before the flux caUbration were around ACaT = 
0.4 A. This value is higher than typical uncertainties due to 
random errors and any other systematic effects (as we will 
see in Section , and stresses the importance of an accurate 
flux calibration before performing any meaningful compar- 
ison between evolutionary synthesis model predictions and 
measured spectra. 



4.4-4 Sky subtraction 

One of the major problems to face in the measurement of 
spectral features in the near-IR is the presence of both, emis- 
sion lines produced by the OH radical (e.g. Rousselot et al. 
2000), and telluric absorptions arising from water vapour 
and other molecules (e.g. Stevenson 1994, Chmielewski 
2000). A proper removal of these effects heavily relies on 
the absence of systematic errors in the data reduction. For 
instance, biases in wavelength calibration (even small sam- 
pling aliasing) and in flux calibration (improper flatfleld- 
ing, residual fringing, and variable slit width along the slit, 
among others) seriously compromise an appropriate sky sub- 
traction. These problems can be specially severe when ob- 
serving faint sources, where the systematic deviations pro- 
duced by an improper sky subtraction likely exceed random 
errors. Since the actual relevance of the mentioned uncer- 
tainties is a function of many variables, including the radial 
velocity of the targets, a case by case analysis must be per- 
formed to cope with them. Note, however, that the avail- 
ability of error frames (in which pixels at the position of 
sky emission lines have larger errors) together with the use 
of error weighted least-squares fits to determine the pseudo- 
continuum (see Section [4.3| ) alleviate the effects of sky resid- 
uals in the measurement of generic indices. 
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Table 8. Calibrations between different calcium index systems. CaT: New calcium index defined in this paper. CaT(j): Calcium index 
in the system by i measured in our spectra. CaT(i)j: Calcium index in the system by i measured in the spectra by j. Codes: JAJ (Jones 
et al. 1984), A&Z (ArmandrofT & Zinn 1988), DTT (Dfaz et al. 1989), ZHO (Zhou 1991), D&H (Delisle & Hardy 1992) and ITD (Idiart, 
Thevenin & de Freitas Pacheco 1997). cr: unbiased standard deviation of the fit. N: Number of stars. Range: Effective temperature region 
where the calibration was obtained. 





Calibrations 




N 


TcS range (K) 


CaT 


= 2.064 H 


- 0.738 CaT(JAJ) 


0.86 


627 


3600- 


38400 


CaT 


=-0.155 H 


- 1.183 CaT(A&Z) 


0.34 


582 


3600- 


-7500 


CaT 


= 1.450 H 


- 1.000 CaT(DTT) 


0.39 


486 


4000- 


-6300 


CaT 


= 0.002 H 


- 1.345 CaT(ZHO) 


0.35 


591 


3000- 


-7500 


CaT 


=-0.756 H 


- 1.107 CaT(D&H) 


0.36 


559 


3600- 


-7000 


CaT(JAJ) 


=-3.322 H 


- 1.524 CaT(JAJ)jAj 


1.15 


21 


3600- 


-6500 


CaT(A&Z) 


= 0.000 H 


- 1.000 CaT(A&Z)iTD 


0.23 


41 


4350- 


-6300 


CaT(DTT) 


=-0.917 H 


- 1.103 CaT(DTT)DTT 


0.55 


102 


3425- 


-6800 


CaT(ZHO) 


= 1.746 H 


- 0.616 CaT(ZHO)zHO 


0.66 


36 


3350- 


-6250 



4.5 Conversions between different calcium index 
systems 

In this section we give a set of calibrations to make conver- 
sions between the different systems of calcium indices. These 
will be useful for readers interested in using our system or 
in transforming their own calcium index data into the new 
indices. 

First, we have compared the measurements of the previ- 
ous calcium index definitions over the 706 library stars with 
the corresponding CaT values. The calibrations were com- 
puted by deriving an error weighted least-squares fit to a 
straight line. If the slope of such fit was not significantly dif- 
ferent from one, we fitted a constant offset. If this term was 
not significant, we kept the one-to-one relation. Due to the 
locations of the sidebands of some previous definitions (see 
Section 4.2), hot and cold stars may depart from the general 



trend of the rest of the library stars. When this occurs, and 
in order to ensure the quality of the fit, we have restricted 
the range of the calibration excluding from the fit those stars 
with extreme temperatures. In particular, hot stars have not 
been used for the comparison with A&Z, DTT, ZHO and 
D&H, where the indices were underestimated due to low 
continuum levels. Similarly, the high continuum levels de- 
rived for the cold stars in the systems by JAJ and DTT lead 
to unreliable indices. Moreover, stars within the T^s range 
but deviating more than Scr from the fitted relation were 
also rejected. Most of them are cluster stars with low S/N 
ratios which present high uncertainties in their Ca measure- 
ments. This becomes specially significant in the comparison 
with JAJ, whose index degrades rapidly as the S/N ratio de- 
creases, producing the large dispersion in the fitted relation. 
Figs. ^ to ^ show the calibrations explained above, includ- 
ing the fitted and rejected points and the final fit. CaT(i) 
refers to the index defined by the reference i measured on 
our spectra. 

Further calibrations are also presented for those papers 
which measure a calcium index over a sample of stars in 
common with our stellar library. Following the same proce- 
dure to compute the fit, the published index values for the 
stars in common were compared with the ones obtained by 
measuring their index over the same stars in our stellar li- 
brary. It must be noted that, since the measured index is 
the same in these comparisons, systematic differences be- 



tween the two sets of measurements are mainly due to flux 
calibration effects (some authors did not perform any ffux 
calibration) and differences in spectral resolution (although, 
given the previous resolutions, see Table ^, this is not the 
dominant effect). The importance of such effects are shown 
in Figs. ^ to |^, where CaT(j)j now refers to the index de- 
fined by i and measured by j. 

In Table ^ we present the derived calibrations together 
with other details of the fits. These fits can be useful for 
those readers who want to convert their data of previous Ca 
indices into the new CaT system. The upper calibrations 
in Table |^ should only be applied once the spectra are on 
the same spectrophotometric system (i.e. equally fiux cali- 
brated) and spectral resolution as the sample of this work. 
If that is not the case, calibrations like the lower ones in 
Table |^ should be applied before. 

Readers interested in using our Ca system to make com- 
parison of observed line-strengths with the fitting-function 
predictions that will be presented in the upcoming papers 
of this series, are highly recommend to convert their obser- 
vations to this system by: 

(i) fiux calibrating the spectra, 

(ii) observing a significant sample of stars in common 
with our stellar library, 

(iii) converting to o ur spe ctral resolution using the rela- 
tions given in Section 1.4.2| and those that will appear in 
Paper IV for the integrated spectra of galaxies. 



5 COMPUTATION OF THE INDEX ERRORS 

One of the principal objectives of this series is to obtain 
a functional representation of the behaviour of the calcium 
triplet as a function of the atmospheric stellar parameters 
(see Paper III). In order to ensure the precision of the fi- 
nal fitting functions, an accurate treatment of the errors is 
needed. To do this, two kind of error sources must be consid- 
ered: the uncertainties in the atmospheric parameters (which 
are thoroughly discussed in Paper II) and those associated 
with the index measurements, which are the subject of this 
section. 

Following the error analysis described in Gorgas et al. 
(1999), we divide the error sources into random errors and 
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CaT(A&Z)„„ CaT(DTT)nTT CaT(ZHO), 



Figure 6. Comparison between different systems. Diagrams (a) to (e) compare the CaT index with the previous indices measured over 
the 706 stars of the stellar library. CaT(i) means the index in the system by i measured over our spectra. Diagrams (f) to (i) compare 
measurements of the same index in our and other spectra for a subsample of stars in common. CaT(i)j means the index in the system 
by i as measured by j. Label codes are JAJ (Jones et al. 1984), A&Z (Armandroff & Zinn 1988), DTT (Di'az et al. 1989), ZHO (Zhou 
1991), D&H (Delisle & Hardy 1992) and ITD (Idiart, Thevenin & de Freitas Pacheco 1997). Symbol types, indicating different ranges of 
effective temperature, are the same as in Fig. M. The dashed line shows the one-to-one relation. Grey symbols are stars deviating more 
than 3o" from the fitted relation, whereas open symbols refer to those stars with effective temperatures outside the fitted range. The solid 
line marks the most significant fit to the symbols in black. In all the cases, indices are measured in angstroms 



systematic effects. In brief, the accuracy of the derived ran- 
dom errors will be checked by the comparison of stars with 
repeated observations within the same run, whereas the sys- 
tematic errors will be mainly derived by comparisons of stars 
in common between different runs. 



5.1 Random errors 



We consider three principal sources of random errors: photon 
statistics, flux calibration and the combined effect of wave- 
length calibration and radial velocity uncertainties. There- 
fore, expected random errors for the i"^ star can be com- 
puted through the following quadratic addition 
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^ [-^a] expected, z " 

photon, i + a [Xa] + a [Xa] wavelength ,i 7 



(8) 



where Xa refers to any of the indices. 

However, unknown sources of random errors may be 
present in the data. Following the method described in 
Gonzalez (1993), we checked whether the standard devia- 
tion of index measurements of stars with multiple observa- 
tions within the same run was significantly larger than the 
expected error, o"[Xa]oxpcctcd, for the whole run (using a F- 
test of variances with a significance level a = 0.1). When 
this occurs, a residual random error cr[Xa]rcsiduai(i) was de- 
rived and introduced in the final random error of each star, 
o"[Xa] random, i, through a multiplicative factor / 



^ [Xa]random,z ^ [Xa] expected, z 7 

which can be computed as 



[Xa] 



residual{l) 



nXa 



(9) 



(10) 



expected 

(see Appendix B for a detailed justification of the above 
equations) . The multiplicative factor for each run was com- 
puted by using the index CaT. ft was only statistically signif- 
icant for some runs, being always / < 1.7. In particular, run 
2 showed an exceptional agreement between the expected 
and measured CaT errors and no correction was applied. 

In the following, a detailed description of the mentioned 
sources of random errors is given: 

(i) Photon statistics and read-out noise. With the aim 
of tracing the propagation of photon statistics and read- 
out noise, we followed a parallel reduction of data and error 
frames with the reduction package REDmE , which creates 
error frames at the beginning of the reduction procedure 
and translates into them, by following the law of combina- 
tion of errors, all the manipulations performed over the data 
frames. In this way, the most problematic reduction steps 
(flat-field and distortion corrections, wavelength calibration, 
sky subtraction, etc) are taken into account and, finally, each 
data spectrum has its corresponding error spectrum which 
can be used to derive reliable photon errors in the index, 
o"[Xa] photon. Typical errors for the whole sample of stars due 
to photon noise are o" [CaT*] photon = 0.158, cr[CaT]photon = 
0.222 and cr[PaT]photon = 0.140. A detailed description of 
the estimation of random errors in the measurement of clas- 
sical line-strength indices is shown in CGCG, whereas an 
analogous treatment for the new generic indices defined in 
this paper is given in Appendix A2. More details concerning 
the sensitivity of different indices to the signal-to-noise ratio 
have already been explained in Section [l.^ . 

(ii) Flux calibration. During each observing night, we 
observed a sample of spectrophotometric standard stars 
(from Oke 1990) and derived a flux calibration curve for 
each one of them. Within the same run, the indices were 
measured using a final flux calibration curve which was com- 
puted by averaging all the individual ones. In addition, we 
used the individual curves to estimate the random error in 
the flux calibration as the r.m.s scatter among the differ- 
ent index values obtained with each one of them. Such an 
error depends on the observing run, but typical values in- 
troduced by this effect are cr[CaT*]flux = 0.033, (T[CaT]flux 
= 0.059 and cr[PaT]flux = 0.057, being almost negligible in 
comparison with the photon errors. 



(Hi) Wavelength calibration and radial velocity. The 
combined effect of wavelength calibration and radial veloc- 
ity correction is another random error source. Radial ve- 
locities for Coma, Hyades and field stars were taken from 
the Bibliography of Stellar Radial Velocities (Abt & Biggs 
1972) and from the Hipparcos Input Catalogue (Turon et al. 
1992), the later giving mean probable errors of ~ 5 km s~^. 
The adopted values for the cluster stars were the averaged 
cluster radial velocities: M3, M5, MIO, M13, M71, M92 and 
NGC6171 were taken from Hesser, Shawl & Meyer (1986), 
NGC188 from Friel (1989), and M67 and NGC7789 from 
Friel & Janes (1993). Typical radial velocity errors for the 
cluster stars are ~ 15 km s"^ (~ 0.43 A at A 8600 A). 
Once the wavelength calibration procedure was applied to 
the whole star sample, we still found small shifts (a typical 
value is 0.32 pixels - or 9.5 km s~^ - for the JKT runs) be- 
tween spectra within the same run, and also systematic off- 
sets (an averaged value is 0.51 pixels) between different runs. 
The deviations within each run were probably due to uncer- 
tainties in the published radial velocities and random loca- 
tions of the star across the spectrograph slit, along the spec- 
tral direction. In order to correct for these effects, we cross- 
correlated each library star with an accurately calibrated 
reference spectrum (a hot or cool library star depending on 
the spectral type of the problem spectrum). The derived 
shifts were applied to the corresponding spectra obtaining 
thus a highly homogeneous wavelength scale for the whole 
stellar library. In spite of that, the indices were measured as- 
suming a radial velocity error of 5 km s~^ which translates 
into rather small errors, typically cr [CaT*] wavelength = 0.003, 

(T [CaT] wavelength = 0.015 and Cr[PaT] wavelength = 0.003. 

To summarize, photon noise is the main source of ran- 
dom error in the measurement of the present calcium in- 
dices. Note, however, that systemat ic effe cts such as an un- 
proper flux calibraton (see Section 4.4.3) can be the domi- 
nant source of uncertainty. 



5.2 Systematic effects 

The main sources of systematic effects in the measurement 
of spectral indices are spectral resolution, sky subtraction 
and flux calibration. 

(i) Spectral resolution. As it is shown in Section 



1.4, the 



new indices are quite insensitive to broadening within an 
ample range of FWHM values. However, in order to obtain 
an homogeneous spectral library for stellar population syn- 
thesis models, we broadened the spectra in runs 1, 3 and 4 
to the spectral resolution of run 2 (1.50 A). Only the cluster 
stars corresponding to runs 5 and 6 (12 and 15 stars, respec- 
tively) were not broadened. These stars keep the spectral 
resolution given in Table ^. 

(ii) Sky subtraction. The exposure times of most of the 
library stars are short enough to neglect the effect of an un- 
certain determination of the sky level. However, some faint 
cluster stars needed up to 1800 seconds of exposing time, 
and many of those in globular clusters were observed with 
the unavoidable presence of cluster members inside the spec- 
trograph slit which complicated the determination of the sky 
regions. Although there is not a simple recipe to compute 
the relevance of this systematic effect, we are confident that 
the sky subtraction is roughly correct since we do not detect 
unremoved sky lines. In any case, we have implicity quanti- 
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CaT(runs 2 & 3) CaT(runs 1 & 3) CaT(runs 1 & 2) 

Figure 7. Comparison of CaT measurements of stars with observations in the JKT runs (1, 2 and 3). The solid hne shows the error- 
weighted derived offset, whereas the dashed Une (almost hidden by the former fit) displays the one-to-one relation. Error bars for each 
point are also plotted. 



fied the sky subtraction uncertainties, as well as that of all 
the remaining untreated systematic effects, by comparing 
the amplitude of the final random errors of the same stars 
in different runs. T his procedure is explained in greater de- 
tail below (Section 5.3). 

(in) Flux calibration. Systematic errors are expected 
to be present due to possible differences between the spec- 
trophotometric system of each run. In order to guarantee 
that all the spectra are in an homogeneous system, we se- 
lected run 1 as our spectrophotometric reference system 
(since it contains a large number of stars in common with 
other runs) and derived a re-calibration curve for each other 
run by averaging the ones obtained from stars in common 
with run 1. Even so, it is not completely sure that our data 
are in the true spectrophotometric system, so we encourage 
readers interested in using our results to include in their 
observing plan stars in common with our stellar sample in 
order to ensure a proper correction of the data. 



5.3 Final errors 

In order to check that the systematic effects between runs 
have been fully corrected, we have compared the Ca indices 
for the stars in common in the three runs at the JKT apply- 
ing a t test (see Fig. ^). It is apparent that the agreement 
between runs is satisfactory. 

A new test to improve the random errors derived in Sec- 
tion 5.1 is to compare the index measurements of stars in 



common between different runs. Considering again run 1 as 
our reference system, the rest of runs were compared with 
it to prove the reliability of the previously derived random 
errors. If the r.m.s. from the index comparison agreed with 
the expected error obtained from their cr[CaT]random,i values 
(using again an _F-test of variances), the run being consid- 
ered was added to the reference system enlarging the fully 
calibrated set for the following comparison. If this was not 
the case, we proceeded as in Section 5.1 and a second resid- 
ual error cr[CaT]residuai{2) 



was derived and introduced to the 



total random error of each star 

Cr[CaT]i = g Cr[CaT]random,i, 

where g is now computed as 



(11) 



1 + 



0-2 [CaT] 



rcsidual{2) 



(T2[CaT] 

expected 



(12) 



Each time the problem run required a g factor, the whole 
comparison procedure was repeated. Fortunately, that was 
only the case of run 6 {g = 1.65) which converged in one 
iteration, and no additional corrections were needed for the 
other runs. 



6 FINAL INDICES AND THE SPECTRAL 
DATABASE 

After converting the spectra to the same resolution and 
wavelength scale, final spectra for the repeated stars were 
obtained through an error-weighted sum of the individual 
spectra. An electronic database containing the final spectra 
of the 706 stars of the stellar library is available at the URL 
addresses; 



http : //www . ucm . es/inf o/Astrof /ellipt/CATRIPLET . html 

and 



http : // www . nott ingham . ac . uk/ "ppzrf p/ CATRIPLET . html 



To summarize, the spectral library spans the range from 
8348.85 to 9019.50 A with a spectral resolution of 1.50 A 
(FWHM) although a few cluster stars from runs 5 and 6 keep 
the original FWHM given in Table |^. The above web page 
also provides an electronic table listing full information for 
each star. It includes the indices CaT*, CaT and PaT mea- 
sured over the final spectra as well as their corresponding 
errors. The Henry Draper Catalogue number, other names 
(mainly HR and BD numbers), coordinates (R.A. and Dec), 
spectral type, luminosity class, apparent magnitude and at- 
mospheric parameters (that will be derived in Paper II) are 
also given. This electronic table is also available at the CDS 



http : //cdsweb.u-strasbg. f r/Cats .html 



7 SUMMARY 

We present a new stellar library in the near-IR spectral re- 
gion covering the range AA 8348-9020 A at a spectral res- 
olution of 1.5 A. It consists of 706 stars with atmospheric 
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parameters in the range 2750 K < T^s < 38400 K, 0.00 
< logff < 5.12 dex, and -3.45 < [Fe/H] < +0.60 dex (see 
Paper II). The aim of this library is to obtain an accurate 
empirical calibration of the behaviour of the Call triplet 
in individual stars (see the derived fitting functions in Pa- 
per III), as well as to derive reliable predictions of the Ca 
strength and spectral energy distribution of stellar popula- 
tions in a wide range of ages and metallicities (see Paper IV) . 
For this purpose, we have defined a new set of line-strength 
indices, namely CaT*, CaT and PaT, which overcome some 
of the hmitations of previous Ca index definitions in this 
spectral range. In particular, the CaT* index is specially 
suited to remove the contamination from Paschen lines in 
the integrated spectra of galaxies. Also, the different sources 
of random and systematic index errors in these and previous 
indices have be en a nalyzed in detail. 

In Section 4.5 we give some recipes for those readers 
interested in using our system or in converting their own 
calcium index data into the new CaT indices. Final spectra 
for the 706 stars of the stellar library and an electronic table 
listing the index measurements, errors and any other infor- 
mation for each star are available at the web pages given 
in Section ^ At the same location we also include a set of 
subroutines to compute the new indices, together with their 
random errors. 
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APPENDIX A: COMPUTATION OF GENERIC 
INDICES AND ASSOCIATED RANDOM 
ERRORS 

In this appendix we give full details concerning the definition 
and practical measure of a new type of line-strength index, 
which will be referred as generic index. In particular, the 
new indices defined and introduced in this paper, namely 
the CaT, PaT and CaT*, are examples of generic indices. 

First, we briefly summarize the main properties of clas- 
sical indices. Their generalization into generic indices ap- 
pears, in a natural way, when introducing three new require- 
ments in the measure of spectral features. The necessity of 
such a change is readily justified when considering the ad- 
vantages of the new definition. Finally, we face the derivation 
of approximate formulae for the estimation of random errors 
from the mean signal-to-noise ratio per A. 



Al Classical indices 

The aim of the use of line-strength indices is to obtain a 
quantitative measure of a spectral signature. Up to date, 
most authors have employed line-strength indices whose def- 
initions are close to the expression of equivalent width. 



Wx{k) = / [1 - S(A)/C(A)] dA, 



(Al) 



where S'(A) corresponds to the observed spectrum, and C(A) 
is the continuum level at the line location obtained by the 
interpolation of S{X) in the local continuum. Note that, as 
it was pointed out by Geisler (1984) (see also Rich 1988), at 
low spectral resolution a pseudo-continuum is measured in- 
stead of a true continuum. In order to prevent subjective 
determinations of the continuum and absorption regions, 
atomic indices have been established by characterizing each 
index with the help of three wavelength bandpasses: the rel- 
evant spectral feature is covered by the central bandpass, 
whereas the other two bandpasses, located at both sides of 
the central region, are employed to define the continuum 
reference level through a linear interpolation. In particu- 
lar atomic indices are computed as (using the notation em- 
ployed by Gonzalez 1993) 



/a(A) 



[1 - S{\)/C(\)] dA, 



(A2) 



where Acj and Acj are the limits of the central bandpass, 
S'(A) is the observed spectrum, and C(A) is the local pseudo- 
continuum, which is derived by 



C{\) = Sb^^-^ + Sr ^ , where 



S(\) dA 



(Ab2 - AbJ 



Sr 



S{\) dA 



(Ar2 — ^ri ) 



Afc = (Afcj + Af,2)/2, A,. = (Ari + Ar2)/2, 



(A3) 



(A4) 
(A5) 



being Af,j , A^j , A^^ , and Arj the limits of the blue and 
red bandpasses respectively. Following the same philosophy, 
broader spectral lines and molecular-band features are usu- 
ally measured in magnitudes using the so-called molecular 
indices, which are computed as 



7m(mag) = -2.5 log^ 



Ac 



(A6) 



The error estimate in atomic and molecular indices was 
extensively analysed in CGCG (see also Cardiel 1999). These 
authors derived a set of accurate formulae for the compu- 
tation of random errors in line-strength indices. The appli- 
cability of the error formulae depends on the availability of 
reliable error spectra. For this purpose, error frames can be 
created at the beginning of the reduction procedure and pro- 
cessed in parallel with data frames. In addition, it is possible 
to use simple expressions to estimate the absolute index er- 
rors as a function of the mean signal-to-noise ratio per A. 
In particular, the error of the atomic indices can be derived 
by 

^ ^ SN{k) ^ ' 

whereas for the molecular indices the expression is 
2.5 logio e C2 1.086 C2 



ci 



(A8) 



SN<yk) SN{A) SN{A)' 

bemg SN{A) the averaged signal-to-noise ratio measured in 
the pixels included in the three bandpasses, and Ci and C2 
two constants which depend on the index definition. 

A2 Generic indices 

The classical index definition can be generalizated by in- 
cluding several spectral features in a single measurement. In 
particular it is possible to define a generic index which can 
be characterized by an arbitrary number of continuum and 
spectral-feature bandpasses. Apart from this modification, 
we have also considered two additional requirements: the 
possibility of modifying the contribution of each spectral- 
feature bandpass (by defining a multiplicative factor for 
each spectral feature), and the derivation of the pseudo- 
continuum using an error weighted least-squares fit. 

Taking all the mentioned requirements into account, the 
generic index can be defined as 



S{X)/C{\)] dA 



(A9) 



where the subscript "a" indicates that this index is the 
generalization of the classical atomic index (i.e. it is mea- 
sured in A), Nf is the number of spectral- feature band- 
passes, ^{k) is the multiplicative factor associated to each 
of these bandpasses, and Aci(fc) and Xc2{k) are the limits of 
the fe"" spectral-feature bandpass. The pseudo-continuum is 
assumed to be obtained through an error weighted least- 
squares fit to Nc continuum bandpasses. 

In practice, the integrals in Eq. (At) must be replaced 
by summations, i.e. 



iVt f N(k) 

Uk) ^eY.\ E [1 



S{\k,i)/C{\k,^)] 



(AlO) 



where O is the linear dispersion (in A/pixel, assumed to be 
constant along the spectrum), N{k) is the number of pixels 
covering the k*^ spectral-feature bandpass, and A^^^ is the 
central wavelength of the i^^ pixel in the k^^ bandpass. The 
continuum level is derived as (see e.g. Bevington 1969) 
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with 



Ql = — {E3 S4 — E2 E5} , 



(All) 



(A12) 



cy-2 — {El E5 — E2 E4} . 



(A13) 

A = El E3 - E2 E2, (A14) 
where we have defined the following parameters 

JVc M(n) 

(A15) 



El 



n — 1 /I — 1 



E2=2. 



\ ^ \ ^ A„,h 



Ti — 1 /l — 1 



X2 



JVc A/{n) 

n = l h = l 
JVc A/{n) 



n = l h = l 



(A16) 
(A17) 
(A18) 

(A19) 



E4 = 
and 

JVc M{n) 
„ _ \ ^ \ ^ Xn,h S(Xn.h) 

n — 1 Ai — 1 

being M{n) the number of pixels covering the n**" contin- 
uum bandpass, and a-^[S{X„^h)] the variance of the observed 
spectrum at the central wavelength of the h^^ pixel in the 
n*** continuum bandpass. 

It is important to note that since, in general, the band- 
passes limits will not be coincident with the wavelengths 
corresponding to the edges of each pixel, fractions of pixels 
must be considered when performing the summations at the 
borders of such bandpasses. 

Using the same notation, the generalization of the clas- 
sical molecular index is straightforward, namely 

Tm(mag) = 



-2.5 log 



1 - 



(A20) 



where the subscript "m" denotes that this index is measured 
in magnitudes. 

Next, and in order to derive analytical expressions for 
the computation of random errors in the generic indices, we 
follow a similar treatment to that described in CGCG for the 
classical indices. The expected random error in the measure- 
ment of a generic index of the form given in Eq. (AlO) can 
be expressed as 



e 



JV(fc) 



E 1 ^(^) E [5(Afe,,)/C(A, 



(A21) 



where the right hand term is a function of 2 X X^^Ii ^i^) 
variables, (. . . , S{Xk,i), ■ ■ ■ , C{Xt,i), ■ ■ .), which are not all in- 
dependent. In fact, if one assumes that none of the A^f -I- A'c 
spectral-feature and continuum bandpasses overlap, and 
taking into account that C(A) is a linear fit to iS'(A) in TVc 
continuum bandpasses, the covariance terms verify 



cov(S(Afc,i), C(Amj)) =0 Vfc,m,i,j, 
cov(C(Afc,i), C(Amj)) /O Vfc,m,i,j, 

V k,m,i,j 

except for i — j when k = m, 



cov(S'(Afc,,),5(A™,j)) =0 
'fc,m G [1,N{], 



(A22) 



with 



i,j e [i,N{k)]. 



Wit h thes e results in mind, and taking partial derivatives in 
Eq. (A21), it is easy to obtain 



e 

JVf JV(i) 



,2,,, C'(Ai,,) a''[S{Xi,i)] + g^(A;,0 a'[C(Ai,i)] 



i=i 1=1 ^ 

JVf JV(i) JVf JV(m) 

EEE E ■ 

i = l i=l m=l 3 = 1 

j^i when m = ( 



cov(C(Aj,0,C(A™,,)) 



r (A23) 



where (J^[C{Xk,i)] is the variance of the fitted continuum, 
evaluated at the 

i*'' pixel of th e k^^ spectral-feature bandpass. Remember- 



ing that Eq. ( All ) is a function of 3 x '^^1-^ M{n) vari- 
ables, (. . . , A;,r, . . . , S{Xi,r), . . . , o-s(A,^), . . .), of which only 
(. . . , ^(A;,^), . . .) have error, it is possible to write 



JVc M(l) 

a'[CiX,,)]^Y.T. 



j dC{Xk,i) \ _2 



\dS(Xl,r) 



^'[^(Aj,.)], (A24) 



with 

ac(Afc..) 

dS{Xl,r) 



A V'^'[5'(A,,,)] 

Ai,r 



Xl.r 



Xk,i 



^HS{Xl,r)] 



El 



<^'[S(A,,.)] 
1 

" <t2[S(A,,0] 



+ 

(A25) 



The non-null covariance terms in Eq. (A2J) can be derived 
using the definition of covariance 

COv(C(Afe^i), C(Amj)) = 

(C(Afc,0 C(A„,,)) - (C(Afc,0) (C(A™,,)) = 

[(aiai) - (qi>(qi)] + (A26) 

[{0102) - (ai)(a2)l {XkA + Amj) -I- 

[{0202) — (a2)(a2)l Xk,i Xm,j- 



Using the parameters defined in Eqs. (A15)-(A19), and tak- 
ing into account that their variances verify 

a^pi] =a2[E2] ^a^pa] =0, 

a2[E4] = El, and cr^ [E5] = E3, 

is not difficult to show that 



(Qiai) - (ai){ai) = ^[Ei E3 E3 - E2 E2 E3] 



(aia2) - (ai){a2) 



-[E2 E2 E2 — El E2 E3 



(A27) 

(A28) 
(A29) 
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SN (A) 

Figure Al. Comparison of absolute random errors in the 
measure of CaT, PaT and CaT* in the 706 stars of the li- 
brary, obtained with the new formulae and with numerical 
simulations (1000 simulations per star and index were per- 
formed), as a function of the mean signal-to-noise ratio per A 
in the bandpasses of each index. There is an excellent agree- 
ment, within statistical errors, between both methods (the mean 
simulations /f^f-^a] formulae value for the 2118 points is 0.99, 
with a r.m.s. of 0.02). 



{Q2Q2) - {012){012) = ■^Pl ^1 ^3 - Si E2 E2]. 



(A30) 



On the other hand, the expected ran dom error in 
generic indices measured in magnitudes, Eq. (A2C), can be 
computed with 



2.5 logipe 



10 



-0.4 Jm 



(A31) 



We have verified the reliability of the new set of formu- 
lae by comparing the measured CaT, PaT and CaT* random 
errors in the 706 stars of the library, with the errors de- 
rived with the help of numerical simulations (see Section 4 
in CGCG). This comparison, which is graphically displayed 
in Fig. Al, indicates that agreement is complete. Note, how- 



ever, that the use of analytical formulae instead of numerical 
simulations should be preferred since the later constitutes a 
computer-time demanding method. 

When the previous formulae are applied over spectra 
corresponding to redshifted objects, it is important to re- 
member that generic indices of the form given in Eq. ( A9), 
and their errors, must be corrected to rest frame values by 



1 + z 



(A32) 



where Ta\^ and Xa|2=o ^''^ redshifted and rest frame 
indices, cr[Xa]z and a[Xa]z=o are their corresponding errors, 
and z is t he redshift. Generic indices of the form given by 
Eq. (A20), in magnitudes, and their errors, do not require 
correction by redshift. 

It is important to note that an appropiate treatment 
of covariance terms is essential in order to guarantee the 
accuracy of the error estimate. Neglecting such terms can 
lead to important systematic anomalies in the determination 
of random errors. For illustration, we have studied this effect 
on the generic indices CaT, PaT and CaT* . For this purpose, 
we measured the individual absorption features of CaT and 



PaT using a new set of indices, CaTi, CaT2 and CaTs, and 
PaTi, PaT2 and PaTa, respectively, which were defined as 
generic indices with the same continuum bandpasses that 
the CaT and PaT indices. In this way, the CaT, PaT and 
CaT* are obviously derived as a linear combination of the 
new set, i.e. 



CaT =CaTi + CaT2 + CaT3, 
PaT = PaTi + PaT2 + PaTs, 
CaT* = CaTi + CaT2 + CaTs- 

0.93 (PaTi + PaT2 + PaTs) 



(A33) 



Using Eq. (|A23D for the computation of random errors in 
CaTi and PaTi (with i = 1, 2, 3) in the 706 stars of the li- 
brary, and adding errors quadratically neglecting covariance 
terms in the set CaTi, PaTi, leads to 

a[ELiCaT,] = (^^^^a2[CaT,])'/', 

a[ELiPaT,] = (2ti'^'[PaT,l)'^', 
a[E?^iCaT, - 0.93E?=iPaTi] = 



(A34) 



(Ell '^'[CaT.] + 0.932 ^3^^ a^paT,]) 



1/2 



The application of these approximate formulae, in c ompar- 
ison with the accurate expression given in Eq. (A23), gives 



(7[E?^iCaT,]/cr[CaT] = 0.7865, r.m.s. = 0.0042, 
(j[E?^iPaT,]/(7[PaT] = 0.9280, r.m.s. = 0.0037, 
(j[E?^iCaT, - 0.93E?=iPaTi]/a[CaT*] = 1.1397, 

r.m.s. = 0.0078. 



(A35) 



Sum mari zing, covariant terms, computed as given by 
Eq. (A26), must be retained in order to guarantee that ran- 
dom errors are free from systematic uncertainties. 



A3 Estimation of random errors from 
signal-to-noise ratios 

As we have already mentioned, it is possible to obtain esti- 
mates of the absolute errors of the classical indices as func- 
tions of the signal-t o-n oise ratio, using expressions as simple 
as Eqs. (A7) and (A8). These formulae constitute a useful 
tool for observation planning. Unfortunately, the approxi- 
mations employed to derive such expressions (see details in 
Cardiel 1999) can not be used in the case of generic indices. 
For this reason we have decided to use an empirical approach 
to this prob lem . 

In Fig. A2 (top panels) we represent the absolute er- 
rors in CaT, PaT and CaT* as a function of the average 
signal-to-noise ratio per A, for the 706 stars of the library. 
An apparently tight correlation is shown by the three in- 
dices, as it is expected. The residuals of a least-squares fit 
(middle panels) reveal that the dispersions of errors (r.m.s.) 
from the linear fit are 3.3, 2.6 and 2.3 per cent for the CaT, 
PaT and CaT* indices, respectively. Although this linear be- 
haviour is a good first estimation of random errors for the 
generic indices, taking into account that the error estima- 
tion of classical atomic indices depends on the index value, 
we have investigated whether a relationship of the form 
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SN{A) 



(A36) 



can also be employed with the generic indices. For this pur- 
pose, we have represented in Fig. A3 the product SN{A) x 
o"[Xa] as a function of the index value, for the stellar library. 
It is clear that a good linear relation is exhibited by the 
three measured indices. The result of a least-squares fit to 
the data is 

18.09 - 0.1751 CaT 



cr[CaT(A)] 
cr[PaT(A)] 
cr[CaT*(A)] 



SiV(A) 

14.27- 0.1463 PaT 

SN{k) ' 

16.43 - 0.1052 CaT* 



(A37) 



SN{A) 



Since the fitted Ci coefficients are two orders of magnitude 
larger than C2 , the dependence of generic index errors on the 
index values is not str ong. However, the application of the 
rela tion s given in Eq. (A37) reduces the residuals showed in 
Fig. 



A2 



ci|^ = (l + z) '- Cil^^o, 
C2\^ = {l+Z)-'/^ C2| 



(bottom panels) for the CaT and PaT indices (r.m.s. 
1.3 and 0.4 per cent, respectively), although it maintains 
those for the CaT* (r.m.s. 2.3 per cent). 

The variation of the coefficients ci and C2 with redshift, 
for classical indices, is described by (Cardiel 1999) 

(A38) 

lz=0 ' 

where ci|^ and C2|^ are the coefficients to be employed in 
spectra with redshift z, whereas ci|^^p and 02 1^^0 are the 
corresponding values at zero redshift. We have verified that 
this is also the case for the generic indices CaT, PaT and 
CaT* (for this purpose, we translated the 706 stellar spec- 
tra of the library as a function of fictitious redshifts, while 
maintaining the signal-to-noise ratio per A ). Therefore, the 
inclusion of the redshift effect in Eq. ( [A36| ) leads to 

(A39) 



(1 + ^) 



1/2 ^1 lz=0 



C2 



T 



SN{A) 



where we have made use of the index variation with redshift 



given in Eq. (A32). Note that although the absolute error 
increases with redshift, the contrary is true for the relative 
error, i.e. 



a[Ja 



= {1 + Z) 



-1/2 



all. 



(A40) 



It is i mportant to highlight that the relations given in 
Eq. (A37) should only be used to obtain approximate es- 
timates of the generic index errors, and that they should 
not be considered as a substitution of the accurate formulae 
described in the previous section. In fact the application of 
approximated equations can lead to misleading index errors 
(see example in CGCG). 



APPENDIX B: RESIDUAL RANDOM ERRORS 

When the unbiased standard deviation of index measure- 
ments of stars with multiple observations, CT[Xa]r.m.s., is sig- 



CaT : 

r \ 1 

: \ : 


PaT : 

: \ : 


CaT* : 

^ \ : 
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Figure A2. {top panels) Absolute errors in the generic indices 
CaT, PaT and CaT* , measured in the 706 stars of the library, as a 
function of the signal-to-noise ratio per A. {middle panels) Resid- 
uals of a least-squares linear fit (rejecting data iteratively outside 
the 99.73 per cent confidence level — indicated by the dashed 
lines — ). {bottom panels) Corrected residuals, obtained b y cor - 
rccting the initial residuals using the relations given in Eq. (A.37). 
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Figure A3. R pipir inal estimation of the constant factors ci and 
C2 of the Eq. (A36). The scatter for low CaT* data are due to 



stars with comparable CaT and PaT indices, which translates into 
relatively small CaT* values, but errors corresponding to indices 
with a larger absolute value. The least-squares fits to straight 
linos have been performed rejecting data iteratively outside the 
99.73 per cent confidence level. 
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nificantly larger than the expected typical error for the 
whole run a"[Ta]expectod, we derive a residual random error, 
o"[Xa] residual by the classical quadratic addition of errors 



[-^ajrcsidual — ^ [^ajr.m.s. ^ [^ajcxpcctcd • (-^-^) 



Although a straight procedure would be to add quadrati- 
cally the residual term to the random error of each star, we 
preferred to introduce it through a multiplicative factor in 
order to keep the original relative qualities among the stars, 
getting at the same time reliable error spectra that can be 
used for any index definition. The introduction of such a 
factor is justified since, as we have shown in Appendix A3, 
random errors of the generic indices CaT, PaT and CaT* can 
be derived from analytical functions of the index value and 
the signal-to-noise ratio (see Eq. (A37)). Thus, we consider 
that the effect of increasing the random errors is equivalent 
to enlarging the photon noise Nk of each star to a new value 
Nl by multiplying with a factor / > 1 

Nk=fN'k, and therefore (B2) 
SNt{k) = SNkik) (B3) 



where k refers to the k^^ star. By substituting Eq. (Bi) in 
Eq. (A36), it can be shown that the final random error, 
o"Pa]flnai,fc, cau be computed as 

final, fc expected, fc ■ 

(B4) 

Next, and in order to derive an analytical expression of /, the 



quadratic addition of errors and Eq. (B4) are simultaneously 
considered 

^ [-^ajfinal.fc ~ ^ [-^ajrcsidual ^ [-^a] expected, fc — 
— ./ ^ Pa]expected,fc, (-B5) 

which leads to a initial expression of / 



^'[^-]-"^-' . (B6) 

^ [-^ajexpectcdjfc 



Note that Eq. (B6) is not definitive since it depends on the 
expected random error of each star. An approximated gen- 
eral expression is finally obtained by substituting the indi- 
vidual errors by the expected typical error for the whole run 



